OPTIMIZED SCHWARZ METHODS FOR MAXWELL'S EQUATIONS 

V. DOLEANf M.J. GANDER+ AND L. GERARDO-GIORDA* 

Abstract. Over the last two decades, classical Schwarz methods have been extended to sys- 
tems of hyperbolic partial differential equations, using characteristic transmission conditions, and 
it has been observed that the classical Schwarz method can be convergent even without overlap in 
certain cases. This is in strong contrast to the behavior of classical Schwarz methods applied to 
elliptic problems, for which overlap is essential for convergence. More recently, optimized Schwarz 
methods have been developed for elliptic partial differential equations. These methods use more 
effective transmission conditions between subdomains than the classical Dirichlet conditions, and 
optimized Schwarz methods can be used both with and without overlap for elliptic problems. A 
simple computation shows why the classical Schwarz method applied to both the time harmonic and 
time discretized Maxwell's equations converges without overlap: for a given frequency we obtain the 
same convergence rate as for an optimized Schwarz method for a scalar elliptic equation. Based 
on this insight, we show how to develop an entire new hierarchy of optimized overlapping and non- 
overlapping Schwarz methods for Maxwell's equations with greatly enhanced performance compared 
to the classical Schwarz method. We also derive for each algorithm asymptotic formulas for the 
optimized transmission conditions, which can easily be used in implementations of the algorithms 
for problems with variable coefficients. We illustrate our findings with numerical experiments. 
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1. Introduction. Schwarz algorithms have experienced a second youth over the 
last decades, when distributed computers became more and more powerful and avail- 
able. Fundamental convergence results for the classical Schwarz methods were derived 
for many partial differential equations, and can now be found in several authorita- 
tive reviews, see [3JI32123], an d books, see SO]- The Schwarz methods were 
also extended to systems of partial differential equations, such as the time harmonic 
Maxwell's equations, see [HI [8], or the time discretized Maxwell's equations, see [39] . 
or to linear elasticity [TH1 [20]) but much less is known about the behavior of the 
Schwarz methods applied to hyperbolic systems of equations. This is true in par- 
ticular for the Euler equations, to which the Schwarz algorithm was first applied in 
[321 133] , where classical (characteristic) transmission conditions are used at the in- 
terfaces, or with more general transmission conditions in [7]. The analysis of such 
algorithms applied to systems proved to be very different from the scalar case, see 

pang. 

Over the last decade, a new class of overlapping Schwarz methods was devel- 
oped for scalar partial differential equations, namely the optimized Schwarz methods. 
These methods are based on a classical overlapping domain decomposition, but they 
use more effective transmission conditions than the classical Dirichlet conditions at 
the interfaces between subdomains. New transmission conditions were originally pro- 
posed for three different reasons: first, to obtain Schwarz algorithms that are con- 
vergent without overlap, see [29j for Robin conditions. The second motivation for 
changing the transmission conditions was to obtain a convergent Schwarz method 
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for the Helmholtz equation, where the classical overlapping Schwarz algorithm is 
not convergent. As a remedy, approximate radiation conditions were introduced in 
[TUl [T2] . The third motivation was that the convergence rate of the classical Schwarz 
method is rather slow and too strongly dependent on the size of the overlap. In 
a short note on non- linear problems [27], Hagstrom et al. introduced Robin trans- 
mission conditions between subdomains and suggested nonlocal operators for best 
performance. In [4], these optimal, non-local transmission conditions were developed 
for advection-diffusion problems, with local approximations for small viscosity, and 
low order frequency approximations were proposed in [301 [S] • In [36j , one can find 
low-frequency approximations of absorbing boundary conditions for Euler equations. 
Independently, at the algebraic level, generalized coupling conditions were introduced 
in [5H1 for discrete overlapping Schwarz methods. Optimized transmission con- 
ditions for the best performance of the Schwarz algorithm in a given class of local 
transmission conditions were first introduced for advection diffusion problems in [SB] , 
for the Helmholtz equation in 6, 25 , and for Laplace's equation in [18]. For complete 
results and attainable performance for a symmetric, positive definite problem, see 
[21], and for time dependant problems, see [231 [22]. The purpose of this paper is to 
design and analyze a family of optimized overlapping and non-overlapping Schwarz 
methods for Maxwell's equations, both for the case of time discretized and time har- 
monic problems, and to provide explicit formulas for the optimized parameters in the 
transmission conditions of each algorithm in the family. These formulas can then eas- 
ily be used in implementations for Maxwell's equations with variable coefficients. As 
we will see, one member of this family reduces in the case of no overlap and constant 
coefficients to an algorithm in a curl-curl formulation of Maxwell's equations, pro- 
posed in [Tj based on [5] , which already greatly enhanced the performance compared 
to the classical approaches in [TJl [8J . 

This paper is organized as follows: in Section [21 we present Maxwell's equations 
and a reformulation thereof with characteristic variables used in our analysis. In 
Section [3j we treat the case of time harmonic solutions. We show that the classical 
Schwarz method for Maxwell's equations, which uses characteristic Dirichlet transmis- 
sion conditions between subdomains is convergent even without overlap. Exploiting a 
parallel with an optimized Schwarz method applied to an Helmholtz equation allows 
us to develop an entirely new hierarchy of optimized Schwarz methods for Maxwell's 
equations with greatly enhanced performance, both with and without overlap. Similar 
equivalence has been presented in [13j for the Cauchy-Riemann equations. In Section 
[U we present and analyze the corresponding hierarchy of optimized Schwarz methods 
for time discretizations of Maxwell's equations. We then show in Section [5] numerical 
experiments in two and three spatial dimensions, both for the time harmonic and 
time discretized case, which illustrate the performance of the new optimized Schwarz 
methods for Maxwell's equations. We also include as an application the cooking of a 
chicken in a microwave oven, a problem with variable coefficients. In Section [SI we 
summarize our findings and conclude with an outlook on future research directions. 

2. Maxwell's Equations. The hyperbolic system of Maxwell's equations de- 
scribes the propagation of electromagnetic waves. It is given by 

- e^- + curl H - a£ = J. + curl £ = 0, (2.1) 

at at 

where £ — (£i,£2,£d,) T and Ti = (J~l\,'H2,'Hz) T denote the electric and magnetic 
fields, respectively, e is the electric permittivity, fj, is the magnetic permeability, a is 
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the electric conductivity and J is the applied current density. We assume the applied 
current density to be divergence free, that is div J = 0. Denoting the vector of physical 
unknowns by 



u = (£ 1 ,£2,£3, Hi, H2, H3) , 
Maxwell's equations (|2.ip can be rewritten in the form 

(G + G d t )u + G x d x u + G y d y u + G z d z u = (J; 0), 
where the coefficient matrices are 



(2.2) 
(2.3) 



G 



oh 



3 



Go 



eh 



fJ-h 



G, 



N, 



x,y,z, 



where O3 (resp. h) represent the 3x3 zero (identity) matrix, and the matrices Ni, 
I = x,y,z are given by 



N,, 




1 
0-10 







1 



-1 






N, 



1 
-10 




For any unit vector n = (m, ni-, 713), ||n|| = 1, we can define the characteristic matrix 
of system (|2.3[) by 



C(n) = Go" 1 



"1 



+ ri2 



N„ 



"3 



-JV, 



whose eigenvalues are the characteristic speed of propagation along the direction n. 
A direct calculation shows that the matrix C(n) has real eigenvalues, 



Al,: 



A 



3,4 



0, 



A 



5,6 



with c = — h= being the wave speed. This implies that Maxwell's equations are hyper- 
bolic, since the eigenvalues are real, but not strictly hyperbolic, since the eigenvalues 
are not distinct, see |2]. For the special case of the normal vector n = (1,0,0), which 
we will use extensively later, we obtain 



G(n) 



■-N x 



whose matrix of eigenvectors is given by 



L 












1 








-Z 











z 








z 











-z 








1 














1 











1 


1 











1 






where Z = y^T denotes the impedance. This leads to the characteristic variables 
w = (wi,W2,W3,W4,W5,'We) T = L~ 1 u associated with the direction n, where 



Wl 

W4 



n 3 



£1, 



W 2 
W 5 



2 (z £ 3 



■n 2 ) 
n 3 ) 



w 3 
w 6 



■■Hi 



\{^ l i-H 2 ). 
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In the following, we will denote by w+, Wq and the characteristic variables asso- 
ciated with the negative, zero, and positive eigenvalues respectively, that is 

w- = (wi, w 2 ) T , w = {w 3 ,W4,) T , w + = (w 5 ,w 6 ) T . (2.5) 

Imposing classical or characteristic boundary conditions on a boundary with unit 
outward normal vector n — (1,0,0) means to impose Dirichlet conditions on the 
incoming characteristic variables W-. for a general normal vector n, this is equivalent 
to imposing the impedance condition (see [2]) 

£ 

B n {£,n) := n x - +n x (H x n) = s. (2.6) 

3. Time Harmonic Solutions. Time harmonic solutions of Maxwell's equa- 
tions are complex valued static vector fields E and H such that the dynamic fields 

£(x, t) = Ke(E(x) exp(iwt)), H(x, t) = Ke(H(x) exp(iwi)) 



satisfy Maxwell's equations (|2.1[) . The positive real parameter uj is called the pulsation 
of the harmonic wave. The harmonic solutions E and H satisfy the time-harmonic 
equations 

— iueE + curl H - oE = J, ilu^iH + curl E = 0. (3.1) 

3.1. Classical and Optimized Schwarz Algorithm. We consider now the 
problem (|3 . If) in a bounded domain f2, with either Dirichlet conditions on the tangent 
electric field, or impedance conditions, on dfl, in order to obtain a well posed problem, 
see [31]. In order to explain the classical Schwarz algorithm for Maxwell's equation, 
we decompose the domain into two overlapping subdomains fix and 0,2, as illustrated 
in Figure 13.11 The generalization of the algorithm formulation to the case of many 
subdomains does not present any difficulties. The classical Schwarz algorithm then 
solves for n = 1,2... the subdomain problems 

-iujeE 1 ' 71 + curl H l n - aE x ' n = J in fi x 



iuj/j,H 1>n + curl E hn = in fii 

L<„ .:/•:"•.//"•; - Bn^E^-^H 2 ^- 1 ) on T 12 

-iueE 2 ' 71 + curl ff 2 '™ - aE 2 ' n = J in fl 2 

iujfxH + curl E = in Q 2 

B n2 (E 2 ' n ,H 2 n ) = B n , 2 (E 1 < n ~ 1 , H 1 '™' 1 ) on T 21 , 



(3.2) 
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where T 12 = dili n S7 2 , ^21 — dil 2 H Q\ and Brij, j — 1,2, denotes the impedance 
boundary conditions defined in (|2.6p . On the physical part of the boundary the given 
boundary conditions are imposed. While the choice of transmission conditions Bn i is 
natural in the view of the hyperbolic nature of the problem, we will see in our analysis 
that there are better choices for the performance of the algorithm. This leads to the 
so called optimized Schwarz methods, 



-iujeE^ 1 + curl H l ' n - aE 1 ' 71 = J in 2 

iujpH 1 ' 71 + curl E 1,n = in 2 

(B ni +SiBn 2 )(E 1 ' n ,H 1 ' n ) = {B Ul + ^BnJ^ 2 '™" 1 , if 2 '"" 1 ) on T 12 
-iujeE 2 ' 71 +cm\H 2 - n - aE 2 ' 71 = J in Vl 2 

iuifxH + curl E = in SI2 

(6„ 2 +5 2 Sn 1 )(-E 2 '", J ff 2 ' n ) = (B n ^+S 2 Bn 1 )(E 1 ' n -\H 1 ' n - 1 ) on r 21 , 

(3.3) 

where 5,-, j = 1, 2 is a tangential, possibly pseudo-differential operator we will study 
in great detail in order to obtain various optimized Schwarz methods. 

3.2. Convergence Analysis for the Classical Schwarz Algorithm. We 

now study properties of the classical Schwarz algorithm (|3.2[) . We use Fourier analysis, 
and thus assume that the coefficients are constant, and the domain on which the 
original problem is posed is = R 3 , in which case we need for Maxwell's equations 
the Silvcr-Muller radiation condition 

lim r(H xn- E) = 0, (3.4) 

r — >oo 

where r — \x\, n = x/\x\, in order to obtain well-posed problems, see |31| . The two 
subdomains are now half spaces, 

fii = (0, 00) x K 2 , fl 2 = (-00, L) x K 2 , (3.5) 

the interfaces are Ti 2 = {L} x R 2 and T 2 i = {0} x K 2 , and the overlap is L > 0. We 
denote by k y and k z the Fourier variables corresponding to a transform with respect 
to y and z, respectively, and \k\ 2 = k 2 + k 2 . 

Theorem 3.1. For any given initial guess (E 1 ' ; if 1 ' ) G (L 2 («i)) 6 , (E 2 '°;H 2 '°) e 
(L 2 (£! 2 )) 6 , the classical Schwarz algorithm &3.ty) with overlap L > 0, including the 
non-overlapping case, is for a > convergent in (L 2 (fli)) e x (L 2 (r2 2 )) 6 , and the 
convergence factor for each Fourier mode k is 



p aa (k,Cj,a, Z,L) 



\J\k\ 2 — Co 1 + iCoaZ — iuj _ v /ffc[2 
-\/|fc| 2 — Co 2 + iuoaZ + id) 



-iuaZL 



(3.6) 



where Co — Lu^/eJI, and Z = is the impedance as before. 

Proof. Because of linearity, it suffices to analyze the convergence to the zero 
solution when the right hand side vanishes. Performing a Fourier transform of system 
(|3.1|) in the y and z direction, the first and the fourth equation provide an algebraic 
expression for E± and H\, which is in agreement with the fact that these are the 
characteristic variables associated with the null eigenvalue. Inserting these expressions 
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into the remaining Fourier transformed equations, we obtain the first order system 

k v k z —Cb 2 +k 2 +iu)n<j 



d x 



[E 2 \ 










E 3 












+ 




lj 2 — k 2 — iiJfUT 


\h 3 J 




iuni 


iuf/j, 
iujfi 



iue-\-<J 
uj 2 — k 2 — iuficr 
iuje-\-a 

o 
o 



iuJS-\-cr 
iuj£-\-a 

o 



(E 2 \ 




(°\ 


E 3 







H 2 







\h 3 J 




W 



The eigenvalues of the matrix in (|3. T|) and their corresponding eigenvectors are 



(3.7) 



\TH 
A 1.2 



and 



-\/|fc| 2 — Co 2 + iujficr, V\ 



(iuie+a)X 
-uj -\-k 2 -\-iujfi(T 
{iuje+a)X 
1 





\ 



v 2 = 



(iuie+a) X 



\ 3 f = \/\k\ 2 - lu 2 + iufia, v 3 = 



V 



* 2 — k 2 — iujficr 
(iw£+cr)A 

l 
o 



o 

i / 

(3.8) 



/ k y —ui +iufj.<7 \ 
(iuie+a)X 



V 4 = 



V 



(itoe-Y-cr) A 

o 
l 



(3.9) 

where we set A := \/|fe| 2 — ui 2 + iujia. Because of the radiation condition, the solu- 
tions of system (|3 . T[) in O;, I = 1, 2, are given by 

(El; E\; Hi- Hi) = (a 1 v 1 +a 2 v 2 )e^ x - L \ (E 2 2 ; E 2 3 ; H 2 2 - flf) = {pivz+frv^e-^ , 

(3-10) 

where the coefficients ctj and f3j (j = 1, 2) are uniquely determined by the transmission 
conditions. At the n-th step of the Schwarz algorithm, the coefficients a — (ai,a 2 ) 
and (3 = (/3±, (3 2 ) satisfy the system 



XL a n-l 



(3 n = B^B 2 e 
where the matrices in the iteration are given by 

k 2 -uj 2 +iuj\ + <jZ{\ + iu>) 



-4, 



A 2 = 



kj-oj 2 +iu)\ + aZ(\ + iu) 



ky k z 



kyk z 



-k 2 +Cj 2 +iuj\ + aZ[\ - iuj) 



k 2 +tu 2 + iuj\ + aZ(X - iQ) 



kyk z 



(3.11) 



and where Bi = Ai, I = 1, 2. 

A complete iteration over two steps of the Schwarz algorithm leads then to 



a 



n+l _ ( a-1 



2 -2AL„.n-l 



{A~{ A 2 ) eT a. 



2-2XL ati-1 



/3" +i = (A^ L A 2 ye- 2AL f3 



and we obtain the iteration matrix 

' \k\ 4 +2XaZ(k 2 -k 2 z ) + X 2 



R={A^A 2 ) 2 e 



2-2XL 



(X+iGi) 2 (X+iOi+crZ) 2 

ik v k z XaZ 
(X+iG>) 2 (X+iO>+aZ) 2 



4k y k z XcrZ 
(X+iCj) 2 (X+iCj+aZ) 2 
\k\ l +2XaZ{k 2 z -k 2 )+X 2 a 2 Z 2 
(X+iCj) 2 (X+iCu+crZy 2 



-2AL 



(3.12) 
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Now by the definition of A, we have |fc| 2 
be re- written in factored form, 



A + w — iujaZ, and thus this matrix can 



i? = 



A — iu 
A + iu 



-2\L 



Id 



AXaZ 



(X + iCb) 2 (X + iQ + aZ)' 2 



ky k z 



ky k-^ 
-k 2 



-2XL 



The convergence factor p c i a of the algorithm is given by the square root of the spectral 
radius of the matrix i?, whose eigenvalues are ( ) e~ 2AL and ( A 7'"T CT f ) e~ 2AL . 

Since a > 0, a direct computation shows that the convergence factor is given by the 
first eigenvalue, which leads to (|3.6[) . and when a ^ 0, a straightforward computation 
shows that p c ia(k) < 1 for all Fourier modes k. □ 
If a = 0, the convergence factor becomes 



p c ia{k,Co,0,Z, L) = 



y/^-\k\2- 



y/u 2-\k\2+ u 

-V|fc| 2 -" 2 £ 



for |fc| 2 < a> 2 , 
for Ifcl 2 > a) 2 . 



(3.13) 



In this case, we obtain for |fe| 2 = Co 2 that the convergence factor equals 1, indepen- 
dently of the overlap, which indicates that the algorithm is not convergent for a = 
when used in the iterative form described here. In practice, Schwarz methods are 
however often used as preconditioners for Krylov methods, which can handle such 
isolated problems in the spectrum. We also see from the convergence factor (|3.13p 
that in the case a — the overlap is necessary for the convergence of the evanes- 
cent modes, |fc| 2 > Cj 2 . Without overlap, L — 0, we have p(|fe|) < 1 only for the 
propagative modes, \k\ 2 < J) 2 , and p(|fe|) = 1 when |fc| 2 > a) 2 . 

Very similar observations were made in the analysis of optimized Schwarz methods 
for the Helmholtz equation in [25] . If one applies to the Hclmholtz equation 



(A + Co 2 )u = f 1 



n = 



— T»3 



(3.14) 



with Sommerfeld radiation conditions limr^oo r (^ — iiuuj = and the same two 
subdomain decomposition p.5p the somewhat particular overlapping Schwarz method 
(note the unequal treatment in the transmission conditions) 



(£ 2 + A)u L { 



f in ^i 
u{ on 1 12 



A)u 



i 

2.U 



f 

(d x - 



-l r (3-15) 
on r 2 i, 



then one obtains precisely the same convergence factor (|3.13p . The classical overlap- 
ping Schwarz algorithm with characteristic transmission conditions (|3.2|) for Maxwell's 
equations is thus very much related to the particular overlapping Schwarz method 
(|3.15p for the Helmholtz problem when a = 0. This particular Schwarz method is a 
very simple variant of an optimized Schwarz method, where one has only replaced one 
of the Dirichlet transmission conditions with a better one adapted for low frequen- 
cies. There are much better transmission conditions for Helmholtz problems, as it was 
shown in [25] . These conditions are based on approximations of transparent boundary 
conditions, which we will study in the next subsection for Maxwell's equations. 

3.3. Transparent Boundary Conditions. To design optimized Schwarz meth- 
ods for Maxwell's equations, we derive now transparent boundary conditions for those 
equations, following the approach in [26] . We consider the time harmonic Maxwell's 
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equations (|3.1[) on the domains fii = (—00, L) x R 2 and £1-2 — (0, 00) x R 2 with right 
hand sides J12 compactly supported in f2i 2, together with the boundary conditions 



(w 2 + +S lW 2 _)(0,y,z) = 0, {w 1 _+S 2 w 1 + ){L,y,z) = 0, (y,z)e 



(3.16) 



and with Silver-Miiller condition on their unbounded part, where w\ and w 2 ^ are de- 
fined in (|2.5p . and the operators Si, I — 1,2, are general, pseudo-differential operators 
acting in the y and z directions. 

Theorem 3.2. If the operators Si, I = 1,2 have the Fourier symbol 



1 



AerZ 



1kyk z 



2 ^ z 

k 2 
K v 



XaZ 



(3.17) 



(A + id;) (A + iu + aZ) 
where A = -\/|fc| 2 — uj 2 + iCoaZ , then the solution 0} Maxwell's equations 13.1}) in f2j 



lOTi/j boundary conditions VS. 10^ coincides with the restriction on fl±.2 of the solution 
of Maxwell 's equations \3.1\) on R 3 . 

Proof. We show that the difference e l , i = 1, 2 between the solution of the global 
problem and the solution of the restricted problem vanishes. We consider the case of 
the second domain, similar computations can be carried out for the first one. This 
difference satisfies in 0,2 the homogeneous counterpart of (|3.1|) with homogeneous 
boundary conditions (13.16|) . and we obtain after a Fourier transform in y and z 



(aiVi + a 2 v 2 )e x + (a 3 v 3 + a^v^e 



-Xx 



where the vectors Vj, j = 1,..,4, are defined in (|3 . 8[) and (|3.9p . The Silver-Miiller 
radiation condition implies that a± = 0:2 = 0. Using now the boundary condition 
(|3.16p at (0, y, z), we obtain that the coefficients ctj, j — 3,4, satisfy the system of 
equations 



0, 



(A 1 +S 1 A 2 ) " 3 

where A\ and A2 arc defined by (|3.11j) . A direct computation 

1 + iibX 



kyk z 



iuj\ 



I/ 2 



ky k z 





«3 




' " 




tt4 








which implies a 3 = 0:4 = 0. Thus e = 0, which concludes the proof. □ 

Remark 1. As in the case of the Cauchy-Riemann equations, see '131, the sym- 
bols in \3. 17\ ) can be written in several, mathematically equivalent forms, 



nsi) = 



{Miu)){\+iuj+aZ) 

1 \—iuj — oZ 71 ,T 



M = 



X—iuj 



M 



\k\ 2 +XaZ A +' w 

(\-iu)(\-iu> - aZ)M~ l , 



where the matrix M is given by 



M = 



XaZ 



2,kyk z 



2kyk z 



k 2 



XaZ 



This motivates different approximations of the transparent conditions in the context 
of optimized Schwarz methods. In the case a = the first form contains a local 
and a non-local term, since multiplication with the matrix M corresponds to second 
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order derivatives in y and z, which are local operations, whereas the term containing 
the square-root of \k\ 2 represents a non-local operation. The last form contains two 
non-local operations, since the inversion of the matrix M corresponds to an integra- 
tion. This integration can however be passed to the other side of the transmission 
conditions by multiplication with the matrix M from the right. The second form 
contains two non-local terms and a local one. We propose in the next section several 
approximations based on these different forms, and analyze the performance of the 
associated optimized Schwarz algorithms. 

3.4. Optimized Schwarz Algorithms for Maxwell's Equations. The trans- 
parent operators Si, I = 1,2, introduced in Subsection 13.31 are important in the de- 
velopment of optimized Schwarz methods. When used in algorithm (|3.3[) . they lead 
to the best possible performance of the method, as we will show in Remark [2] The 
transparent operators are however non-local operators, and hence difficult to use in 
practice. In optimized Schwarz methods, they are therefore approximated to obtain 
practical methods. If one is willing to use second order transmission conditions, then 
the only parts of the symbols in (|3.17p that need to be approximated are the terms 
A = -\/|fc| 2 — u> 2 + iujaZ, because the entries of the matrices are polynomials in the 
Fourier variables, which correspond to derivatives in the y and z direction. 

Theorem 3.3. For the optimized Schwarz algorithm \3. Sty with the two subdomain 
decomposition i3.5\) . we obtain for a = the following results: 
1. If the operators Si and S2 have the Fourier symbol 



ai := T(Si) = 7; 
then the convergence factor is 



h 2 — h 2 

Ky K z 

— Oh h h 2 — h 2 



2k y k z 



li e C(k z ,k y ), 1=1,2, (3.18) 



p = 



(Vife| 2 -s 2 +w) 2 l-TicVifei 2 -" 2 -^) 2 i-72(Vifci 2 -^ 2 -^) 2 

2. If the operators S± and S2 have the Fourier symbol 



(3.19) 



ai := F(Si) = St 
then the convergence factor is 



h 2 — h 2 — Oh h 

hjy n, z Ahjyl^z 

— Oh h h 2 — h 2 

L*f\iyhjg hj^ hjy 



, 7; € C{k z ,k y ), 1 = 1,2, (3.20) 



= (V|fc| 2 -^ 2 +»^) 2 ^-(Vjfcj 2 ^ 2 -^) 2 s 2 -(^\k\*^-icb) 2 -2 X /\k\*-c,n 

(^/\k\ 2 ~^y 2 -iO) 2 5i-(\/|fe|2-<Zi 2 +ii) 2 S 2 -(y/\k\ 2 -0 2 +iCj) 2 

(3.21) 

3. If the operator S\ has the Fourier symbol 13.1 Sty and S2 has the Fourier symbol 
h3.20\) . then the convergence factor is 



1-71 (^Ife!^^^) 2 S 2 -(y/\k\ 2 z u ! 2 r lCj) 2 c -2^/]kT 2 
l- 7l (V|fe| 2 -" 2 -^) 2 6 2 -(^\k\ 2 -u 2 +lQ) 2 



-u 2 L 



1/2 



(3.22) 



Proof. The convergence results are again based on Fourier analysis, as in Section 



13.11 At the n-th step of the Schwarz algorithm, the coefficients a" 
P = {/3i,[3 2 ) in (|3~T0f satisfy 



l ) and 



AT A->e~ 



f3" 



BT 1 ^ e 



XL, 



(3.23) 
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where A = \/|fc| 2 — uj 2 , and the matrices Ai and Bi, I = 1,2, are given by 

A 1 =A l +a 1 A 2 , A 2 =A 2 +a 1 A 1 , B 1 =A 1 +a 2 A 2 , B 2 =A 2 +a 2 A 1 , 

with Ai, I — 1,2, defined in (|3. 1 1|) . A complete double iteration of the Schwarz 
algorithm leads therefore to 

a n+1 = Al 1 A 2 B^ 1 B 2 e- 2XL ct n - 1 , (3 n+1 = B^ 1 B 2 A\ X A 2 e- 2XL . 

Notice that the matrices A\ and A 2 verify the properties 

MA\ = -(A + iu>) 2 A 2 , MA 2 = -{\-iu) 2 A ll (3.24) 

which are essential in all three cases: 

1. We obtain, using ([3T2"4"[) . 

Ax = A 1 + 11 MA 2 = (1 - (A-^) 2 7i)Ai, 

A 2 = A2+J1MA1 = (1- (\ + iQ) 2 ~ fl )A 2 , 

Bi = A1+J2MA2 = (1- (A-iw) 2 7 2 )Ai, 

B 2 = A2+T2MA1 = (1 - {X + iCj) 2 j 2 )A 2 , 

and therefore the iteration matrix becomes 

Rl = A^A 2 B^B 2 = ^^S^^ ^ 1 ^ 2 ^ 

and by using the spectral radius of (|3.12p the result follows. 

2. We get 

Ai = A 1 +8 1 M- 1 A 2 = (l- 
A 2 = A 2 +5 1 M- 1 A 1 = M- 
Bi = A 1 + 82M- 1 A 2 = (l- 
B 2 = A 2 +6 2 M- 1 A 1 = (1- 

and thus the iteration matrix is 

R 2 = A^A 2 B^B 2 = (A±f) 4 ( ^§^^l (AT l A2re-^, 

and we use again the spectral radius of (|3.12|) to conclude. 

3. The conclusion follows as in the first two cases. 



Si > 


\M 


(X+iQ) 2 j 


Si ' 


\A 2 


(X-iCj) 2 . 


S 2 ' 


\M 


(X+iQ) 2 1 


S 2 ^ 


)A 2 


(X-iu) 2 j 



Remark 2. From kS.19\) , we see that the choice 71 = 72 = l/(^/\k\ 2 — Co 2 + iuj) 2 
is optimal, since then p=0, for all frequencies k. With this choice 0/71 and 72, the 
matrices A 2 and B 2 actually vanish. 

3.5. A hierarchy of optimized transmission conditions. We present now 
several particular choices of the remaining parameters in the transmission operators 
Si in Theorem 13 .31 To facilitate the use of our results in domain decomposition codes, 
we return to the initial notation using the physical parameters u>, e and fi. 
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Case 1: taking jx = 72 = in (|3.18p . which amounts to enforce the classical charac- 
teristic Dirichlet transmission conditions, the convergence factor is 



pi(w,e,/i,L, \k\) = 



-UJ 2 €fiL 



\/\k\ 2 — uj 2 ep + iiv^ep I 



In the non-overlapping case, L = 0, this choice ensures convergence only for 
propagative modes, and corresponds to the Taylor transmission conditions of 
order zero proposed in the seminal paper [llj for the Helmholtz equation. 



Case 2: taking = 72 



1 s — iujsii 



in (|3~T5|) or ^ 



in (SB) 



|fc|2 s+iu^/WJl «— = f 11 ~ \k\ 2 -2u> 2 eti+2iui^/WJIs 

and 82 = |fc| 2 — 2cu 2 ep — 2icv«Jeps in (|3 . 20[) with s € C, the convergence 
factor is 



p 2 (w,s,/i,L, \k\,s) 



2 — Lu 2 e/j, 



y/\k\ 2 - LU 2 £[J, + . 



^-2^/\k\ 2 -u> 2 e^L 



which is for L = identical to the convergence factor obtained for optimized 
non-overlapping Schwarz methods for the Helmholtz equation in |25j . 
Case 3: taking 71 = 72 = 77-7: ^ =rr in (13.181) with s 6 C, the convergence 



factor is 

p 3 (uj,e,p,L, \k\,s) 



-uj^efi—iujy'efi 



y \k\ 2 —uj 2 efj,+iujy/ejl 

< P2(u,e,p,L,\k\,s). 



p 2 (w,s,/i,L, \k\,s) 



Case 4: taking 7; 



1 'i-Wpj / = 1 2in (EH]) or 7l = -w— z in 

(|3~Tg|) and <5 2 = |fc| 2 - 2u?ep - 2iLu^/epZs 2 in (|3T20]) with s; e C, I = 1, 2, the 
convergence factor is 



p±{ijj,e,n,L, |fe|,si,s 2 ) 



-LU z €fl—S 



-u^efi— S2 



V|fc| 2 -t- 2 £M+si Vl fc l 2 -^ 2 ^+^2 



-2-v/|fc| 2 -" 2 eMi 



which is for L = identical to the convergence factor obtained for a two 
sided non-overlapping optimized Schwarz method for the Helmholtz equation 
in [231. 



Case 5: taking 7; 



|fc| 2 -2w 2 £^+2iw v /£7Zs ! 



in (|3. 18|) with s/ € C, I = 1, 2, the conver- 



gence factor is 

p 5 (iv,e,p,L, \k\,si,s 2 ) 



VW 2 



- uj * e (i — i uj sje^i 



p 4 (u,e,ii,L, \k\,s 1 ,s 2 ) 



V \ k\ 2 —uj 2 efi-\-iuj^sJI 

< p 4 (uJ,e,p,L,\k\,s 1 ,S2)- 

Except for Case 1, all cases use second order transmission conditions, even though we 
use only a zeroth order approximation of the non-local operator -\/|fc| 2 — ui 2 ep. In the 
cases with parameters, the best choice for the parameters is in general the one that 
minimizes the convergence factor for all |fe| € K , where K denotes the set of relevant 
numerical frequencies. One therefore needs to solve the min-max problems 

min max Pj(uj, e, p, L, \k\, s), j — 2, 3, min max pj(w, e, p, L, \k\, S\, s 2 ), j — 4, 5. 
s6C \k\eK ' s 1 ,s 2 ec \k\eK 

(3.25) 
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with overlap, L = h 



without overlap, L = 



Case 



P 



parameters 



P 



parameters 



1 - ^k\~Co 2 h 
l-2Cj/ii 

1 - 2tC^hi 
l-2i{k\~Cj 2 )Toh^ 



none 



2-h\ 



1 - 



1 

_ l 

V2Ci 



Vh 



V2(k%-cj 2 )i ^ 



Vc 



— 5 T 



-hi 



2" 



(fc 2 -ii 2 )s 

1 

CI 



Pi 



P2 



_ Vcc? 

' ~ V2~Vh 

Vc(k 2 + -ui 2 )i 

\/2~Vh 

3 ! 



3 

h* „ 



(fci-ii 2 )5 -ci 

Pi = "~T 



Table 3.1 

Asymptotic convergence factor and optimal choice of the parameters in the transmission con- 
ditions for the five variants of the optimized Schwarz method applied to Maxwell's equations, when 
the mesh parameter h is small, and the maximum numerical frequency is estimated by fc max = 4- ■ 
Here u) = LU^eJl and Cq = min (fc^_ — uj 2 , u) 2 — k^_) . 



We can choose K = [(/c m i n ,fc_) U (k + , fc max )] 2 , where fc mi „ denotes the smallest fre- 
quency relevant to the subdomain, and k max = % denotes the largest frequency sup- 
ported by the numerical grid with mesh size h, and k± are parameters to be chosen to 
exclude the resonance frequencies. If for example the domain f2 is a rectilinear conduc- 
tor with homogeneous Dirichlet conditions on the lateral surface, the solution is the 
sum of the transverse electric (TE) and transverse magnetic (TM) fields. If the trans- 
verse section of the conductor is a rectangle with sides of length a and b, the TE and 
TM fields can be expanded in a Fourier series with the harmonics sin(^^) sin( 22 ^), 



where the relevant frequencies are |fe| = ny + fr, m,n € N . The lowest one 



is therefore k m i n = 7ry jpr + 5?) an( i if the mesh size h satisfies h — jf = jj, where 
N and M are the number of grid points in the y and z direction, then the high- 
est frequency would be k max = -^j^-- The parameters k± would correspond to the 



frequencies closest to uiyfejl, i.e. = 7ry + -gj and k + — iry + , where 



r + -gr < cjy/ejl < Try "jr + but such precise estimates are not necessary if 
Krylov acceleration is used, see [251 H3] • 

The complete mathematical analysis of the min-max problems (|3 . 25[) is hard, and 
currently open for L > 0. When L = 0, i.e. no overlap, Case 2 and Case 4 are equiva- 
lent to the corresponding optimized Schwarz method for the Helmholtz equation, for 
which theoretical results are available, see [53]. Here, we use asymptotic analysis and 
an equioscillation principle to solve all the min-max problems in p.25[) asymptotically 
as the mesh size goes to zero, in order to obtain compact formulas for the best param- 
eters to be used in our numerical simulations. This leads to the asymptotic formulas 
for the optimized parameters of the form s = p(l — i) and s; = pi(l — i), I — 1,2, with 
p and pi shown in Table 13.11 These results allow us to compare the performance of 
all the optimized Schwarz methods for Maxwell's equations theoretically: we obtain 
a hierarchy of better and better convergence factors starting with Case 1 and ending 
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with Case 5. In addition, the explicit formulas for the optimized parameters can be 
used in order to easily obtain black-box optimized Schwarz methods for Maxwell's 
equations, which would not be possible otherwise. In Section we will furthermore 
verify these theoretical results numerically. 

4. The Case of Time Discretization. If we do not assume the wave to be 
periodic in time, the time domain also needs to be discretized. We consider a uniform 
time grid with time step At, and use an implicit time integration scheme for the time 
derivative in (|2.1|) of the form 

- + curl [H^Hl) = J, 

l M^_Hl +cml ( £" +1 2 +£" ) =0 , 

where the mean value is introduced to ensure energy conservation, see |14) . With this 
time discretization, we have to solve at each time step the system 

— oE — s^/rjE + curl H = J , f-i^/rjH + curl E = g, (4.1) 

where we have set (E,H) := (E n+1 ,H n+1 ), := J ':= J - ^qeE n + 2aE n - 
curl H n , and g = ^TjfiH n - curl E n . 

4.1. Classical and Optimized Schwarz Algorithm. As in the time harmonic 
case, we consider the problem (|4.1[) in a bounded domain Q, with either Dirichlet 
conditions on the tangent electric field, or impedance conditions, on dQ, in order 
to obtain a well posed problem, see [31]. For the two subdomain decomposition in 
Figure 13.11 the classical Schwarz algorithm would at each time step then perform the 
iteration 

—^/rjeE l,n + curl fT 1,n — aE 1 '" = J 1 in 

1,71 + curl E l ' n = g 1 in fij 

Bn 1 (E hn ,H 1 ' n ) = BndE 2 -"- 1 ,!! 2 '"- 1 ) on T 12 
- ^/veE 2 ' n + curl H 2 ' n ~aE 2 > n = J 2 in fi 2 



(4.2) 



^rjfiH + curl E = g 2 in SI2 

Bn 2 (E 2 ' n ,H 2 ' n ) = Bn^E 1 ' 71 - 1 ,^' 71 - 1 ) on F 



21- 



Theorem 4.1. Let ft = M 3 be decomposed into f2i := (— 00, L) x M 2 and 
tt 2 := (0,+oo) x K 2 , L > 0. Then, for any initial guess (E 1 ' ; H 1 ' ) e (L 2 (^i)) 6 , 
(E 2,a ; H 2 ' ) G (L 2 (Q2)) e , the classical Schwarz algorithm \4- 6 <ty with overlap L > 0, 
including the non- overlapping case, is for a > convergent in (L 2 (fli)) 6 x (L 2 (fl2)) e , 
and the convergence factor is bounded by 



R, a = V^^Z^ e -v^V^V^ < 1, (4.3) 

where f) = ne/i. 

Proof. This result follows like in the time harmonic case, simply replacing iuj by 
The convergence factor in Fourier is 



Pcla{\k\) 
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with overlap, L = h 



without overlap, L = 



Case 


P 


1 


1 - 


2a fji \fh 


2 


1 - 


2 6 7/ 6 tl 3 


3 


1 - 


24778/7,4 


4 


1 - 


4 „ 1 1 

2 5 77 10 /j 5 


5 


1 - 


2 6 7712/7, 6 



parameters 



-2^/t, 
r _ 

4r/I V ft 

Vc 

CI 

C5 



parameters 



P : 
P 



none 

1 1 

- 1 

— v 7 ! 



none 

_ VCfji 



Pi 



n 5 

-T-T'P2 



2~Sh1 



V 5 

— i — r 



2 3 11 3 23 1)6 

Pi = — r-, P2 = —i- 

fr_3 h3 



V 



221 
2 



Table 4.1 

Asymptotic convergence factor and optimal choice 
ditions for the five variants of the optimized Schwarz m 
equations, when the mesh parameter h is small, and the 
by fcmax = % ■ Here fj = r/efi. 



4 1 2 3 

Pi — 1 7 P2 — J 

h_5 hj 



of the parameters in the transmission con- 
ethod applied to the time domain Maxwell's 
maximum numerical frequency is estimated 



and the method thus converges for all Fourier modes. To conclude the proof, it suffices 
to take the maximum of the convergence factor over |fc|. □ 

The preceding theorem shows that the classical Schwarz algorithm with Dirichlet 
transmission conditions applied to the time-discretized Maxwell's equations is conver- 
gent for all frequencies |fe|, and that the overlap is not necessary to ensure convergence. 
The classical Schwarz algorithm corresponds in the case a — to a simple optimized 
Schwarz algorithm for an associated positive definite Hclmholtz equation 

(77 - A)u = /, (4.4) 

and from [21] we know that there are much better transmission conditions for such 
problems. We thus propose at each time step the new algorithm 

-y/rjsE 1 ' 71 + curl H l n - aE l ' n = J 1 in fii 

V^/xff 1 '" + curl E x ' n = g 1 in ^ 

(Bn 1 +«S 1 Bn 3 )(£ 1 ' n ,-W 1 '") = (B Ul + S 1 Bn 2 )(E 2 ' n ~ 1 , H 2 ^ 1 ) on T 12 
-^/TjeE 2 ' n + curl H 2 ' n - aE 2 ' n = J 2 in fi 2 

y/rjp.H + curl E = g 2 in £1 2 

(B n2 +S 2 B ni )(E 2 ' n ,H 2 - n ) = (Bn 2 + S 2 B ni )(E 1 ' n ~ 1 , H 1 '"" 1 ) on T 21 . 

(4.5) 

Now Theorem 13.21 Remark [IJ Theorem 13.31 and all the cases in subsection 13.51 hold 
unchanged for the time discretized case of Maxwell's equations upon replacing iui by 
y/fj, so we do not restate these results here. The nature of the associated min-max 
problems (|3 . 25[) however changes fundamentally, and the optimization parameters are 
now real, s = p EM. and si = pi G R, I = 1, 2. For cases 2 and 4 a complete analysis is 
available, see |21J . Using a lengthy asymptotic analysis again, we complete the results 
for the other cases, and show in Table 14.11 the asymptotically optimal parameters to 
use in the time domain case. Again we obtain an entire hierarchy of optimized Schwarz 
methods, with better and better convergence factors from Case 1 up to Case 5. While 
for the time harmonic equations Case 2 and 3, and Case 4 and 5 were asymptotically 
comparable, here all cases are asymptotically different. It is also interesting to note 
a relationship of the optimized parameters for the time domain case with the one 
for the Cauchy-Riemann equations, see [T3]: Case 2 and 4 are identical, since the 
corresponding convergence rates in the two cases are the same, while for Case 1, 3 
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and 5 there is a small difference in the constants, which is due to the additional low 
frequency term in the Maxwell case. The difference appears to be systematic, the 
convergence factor of the Maxwell case is obtained from the convergence factor of the 
Cauchy-Riemann case by replacing h by 2h, while for the optimized parameters one 
has to multiply by 2 in addition to the replacement of h by 2h. 

5. Numerical Experiments. We discretize the equations using a finite volume 
method on a staggered grid, which leads to the Yee scheme in the interior. For the first 
two test cases we consider the propagation in vacuum with s = \i = 1 and a — 0. We 
first show the two dimensional problem of transverse electric waves, since this allows 
us to compute with finer mesh sizes and thus to illustrate our asymptotic results by 
numerical experiments. We simulate directly the error equations, / = 0, on a uniform 
mesh with mesh parameter h, and we use a random initial guess to ensure that all the 
frequency components are present in the iteration. We then show the full 3d case, first 
for a model problem, and then for the application of heating a chicken in a microwave 
oven. 

5.1. Two-dimensional case. We consider the transverse electric waves prob- 
lem (TE) in the plane (x, y, 0). There is no more dependence on z and the components 
E3, Hi i?2 are identically zero. The problem obtained is formally identical to the 
three-dimensional case (|2.3[) . if u — (E\, E%, #3)*, and the matrix 7V V becomes 

and the matrices G x , G y and G v are 

All the analytical results remain valid, we only need to replace |fe| by \k y \, and the 
corresponding quantities in the optimized parameters for both time-harmonic and 
time-discretized solutions. We solve Maxwell's equations on the unit square £1 = 
(0, l) 2 , decomposed into the two subdomains fii = (0, 0) x (0, 1) and f2 2 = (a, 1) X 
(0, 1), where < a < /3 < 1, and therefore the overlap is L — (3 — a, and we consider 
both decompositions with and without overlap. 

In the time-harmonic case, the frequency Q — 2n is chosen such that the rule of 
thumb of 10 points per wavelength is not violated. Table lOl shows the iteration count 
for all Schwarz algorithms we considered, in the overlapping and non-overlapping 
case. The results are presented in the form its(itcM), where its denotes the iteration 
number for the iterative version of the algorithm and Uqm the iteration number for 
the accelerated version using GMRES. 

In Figure 15.11 we show the results we obtained in a graph, together with the 
expected asymptotics. Both on the left in the overlapping case and on the right in 
the non-overlapping one, the asymptotics agree quite well, except for the classical case 
with overlap, where the algorithm performs better than predicted by the asymptotic 
analysis. In the case of the Cauchy-Riemann equations [13 , it was observed that 
certain discretizations of the hyperbolic system can introduce higher order terms in 
the discretized transmission conditions, which can improve the convergence behavior, 
as we observe it here, an issue that merits further study. 

For the time discretized Maxwell's equations we choose fj = 1. Table I5~2l shows the 
iteration count for all Schwarz algorithms we considered, in the overlapping and non- 
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with overlap, L = h 


without overlap, L = 


= 


h 


1/16 


1/32 


1/64 


1/128 


1/16 


1/32 


1/64 


1/128 


Case 1 


18(17) 


27(21) 


46(27) 


71(33) 


-(48) 


-(73) 


-(100) 


-(138) 


Case 2 


16(13) 


16(14) 


17(15) 


20(17) 


28(22) 


36(26) 


50(34) 


68(40) 


Case 3 


10(12) 


12(13) 


14(14) 


16(17) 


31(20) 


40(23) 


56(25) 


81(28) 


Case 4 


17(13) 


17(14) 


20(16) 


22(18) 


26(20) 


28(24) 


33(28) 


38(30) 


Case 5 


10(12) 


12(13) 


14(15) 


17(18) 


41(24) 


53(26) 


63(30) 


73(32) 



Table 5.1 

Number of iterations in the 2d time harmonic case to attain an error tolerance of = 10 — 6 for 
different transmission conditions and different mesh sizes. 



Asymptotic convergence for the overlapping algorithm: time-harmonic case Asymptotic convergence for the non-overlapping algorithm: time-harmonic case 




10" ? 10 2 



Fig. 5.1. Asymptotics for the overlapping (left) and non- overlapping (right) cases for the time 
harmonic eguations. 

overlapping case. We observe that the classical non-overlapping algorithm converges 
only very slowly, the need of optimized methods is evident here. 

In Figure 15.21 we show the results we obtained in a graph, together with the 
expected asymptotics, and there is very good agreement. 



Asymptotic convergence (or lie one 'I Lipping a .lorirhm: lime-domain case Asymptotic convergence for the non-overlapping algorilhm:time-domain case 




Fig. 5.2. Asymptotics for the overlapping 
discretized equations. 



(left) and non- overlapping (right) cases for the time 
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with over 


lap, L - 


= h 


without overlap, L = 


h 


1/16 


1/32 


1/64 


1/128 


1/16 


1/32 


1/64 


1/128 


Case 1 


17 


24 


33 


45 


280 


559 


1310 


2630 


Case 2 


13 


15 


19 


24 


39 


56 


77 


111 


Case 3 


12 


14 


16 


18 


13 


16 


20 


26 


Case 4 


12 


13 


15 


17 


21 


25 


30 


36 


Case 5 


12 


14 


16 


18 


13 


17 


19 


22 



Table 5.2 

Number of iterations in the 2d time discretized case to attain an error tolerance of 10 — 6 for 
different transmission conditions and different mesh sizes. 





with 


overlap, L = h 


without overlap, L = 


h 


1/8 


1/16 


1/30 


1/8 


1/16 


1/30 


Case 1 


19(13) 


29(17) 


46(22) 


-(93) 


-(140) 


-(202) 


Case 2 


14(12) 


19(14) 


23(16) 


48(29) 


69(36) 


98(48) 


Case 3 


16(12) 


18(14) 


21(16) 


65(35) 


80(42) 


166(55) 


Case 4 


15(13) 


19(15) 


22(17) 


38(28) 


60(33) 


104(39) 


Case 5 


16(13) 


18(14) 


21(16) 


70(36) 


80(42) 


176(55) 



Table 5.3 

Number of iterations in the 3d time harmonic case to attain an error level o/10" 6 for different 
transmission conditions and different mesh sizes. 



5.2. Three-dimensional case. We solve now Maxwell's equations on the unit 
cube f2 = (0, l) 3 . We decompose the domain into two subdomains Vl\ = (0, 0) x (0, l) 2 
and VI2 — (a, 1) x (0, l) 2 , with < a < (3 < 1, and L — (3 — a as before. In the time- 
harmonic case, we chose the frequency u) = 27r/3 to satisfy the rule of thumb of 10 
points per wavelength. Table [5~3l shows the iteration count for all Schwarz algorithms 
we considered, both in the overlapping and non-overlapping case. 

The results for the time discretized Maxwell's equations where rj = 1 are shown 
in Table EH 

5.3. A three-dimensional application: chicken in a micro-wave oven. 

We apply now the previous principles to derive an efficient domain-decomposition 
method based on optimized interface conditions to solve a realistic application: heat- 
ing up a chicken in a micro- wave oven, see Figure [5~31 on the left. The computational 
domain is now given by the heating cavity of a Whirlpool Talent Combi 4 microwave 
oven, £1 = [0, 0.32] x [0, 0.36] x [0, 0.20] meters. We impose metallic boundary condi- 
tions (which means a null tangential electric field) on all faces except on the right of the 
oven, where the components of the electric field are the dominant TE10 mode gener- 
ated by the magnetron on a small rectangle of dimensions 0.08 x 0.04. The electric and 
electromagnetic properties of the media are now non-constant in the computational 
domain: inside the chicken, we have an electric permittivity e = 4.43 • 10~ n Fa |~° d ° and 
the conductivity is a = 3 • 10~ n Sic ™ cn3 ) whereas for the air e = 8, 85 • 10~ 12 Fa ^ d3 and 
a = Slc ™ cns , The magnetic permeability is the same for both, fi = 4ir- 10~ 7 Hc °'" Y , and 
the frequency is given by u> — 2ir ■ 2.45 GHz. 

We decompose the microwave oven into 2x2x2 = 8 subdomains of equal size 
on a grid with mesh size h — 0.005, which allows us to solve this problem on a PC, 
where a direct factorization would not have been possible any more. The real part 
of the magnetic field of the solution is shown in Figure 15.31 on the right, and the 
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with 


overlap, L = h 


without overlap, L = 


h 


1/8 


1/16 


1/30 


1/8 


1/16 


1/30 


Case 1 


14 


18 


25 


246 


467 


859 


Case 2 


13 


18 


22 


46 


65 


87 


Case 3 


12 


15 


17 


47 


59 


73 


Case 4 


14 


17 


19 


48 


57 


66 


Case 5 


12 


14 


16 


46 


53 


60 



Table 5.4 

Number of iterations in the 3d time discretized case to attain an error level of 10 — 6 for different 
transmission conditions and different mesh sizes. 




Fig. 5.3. Chicken in a Whirlpool Talent Combi 4 micro-wave oven on the left, and real part of 
the magnetic field in the cooking cavity while heating the chicken on the right. 

intensity (Euclidian norm) of the electric and magnetic field in the oven are shown 
in Figure [53] in three dimensions. Two-dimensional cross sections of the solution are 
shown in Figure 15.51 where we show in each row on the left the electric and on the 
right the magnetic field intensity. One can see from these computational experiments 
why a turntable is so important in a microwave oven: there are hot spots, where the 
intensity of the standing wave is high in the chicken, and other areas, where there is 
very little heating happening. Only a turntable can lead to an approximately even 
heating of the chicken. 

6. Conclusions. We have shown that for Maxwell's equations, a classical Schwarz 
algorithm using characteristic Dirichlet transmission conditions between subdomains 
has the same convergence behavior as a simple optimized Schwarz method applied to 
the Hclmholtz equation, with a low frequency approximation of the optimal trans- 
mission conditions. This relation allowed us to develop easily an entire hierarchy of 
optimized overlapping and non-overlapping Schwarz methods with better transmis- 
sion conditions than the characteristic ones for Maxwell's equations. We illustrated 
with numerical experiments that the new algorithms converge much more rapidly 
than the classical one, and that such algorithms can be effectively used to compute 
an approximate solution for a large scale application. This latter problem contains a 
positive conductivity, variable coefficients and multiple subdomains, a case which is 
not covered by our current analysis. Nevertheless, the algorithm performs well with 
the coefficients derived from the zero conductivity, constant coefficient case. We are 
currently studying the optimization problem with non-zero conductivity, for which 
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Fig. 5.4. Chicken heating in a microwave oven: electric field intensity on the left, and magnetic 
field intensity on the right. 



the equivalence with the Helmholtz equation does not hold any longer. 

The equivalence between systems and scalar equations has already been instru- 
mental for the development of optimized Schwarz algorithms for the Cauchy-Riemann 
equations, and will almost certainly play an important role for other cases. For exam- 
ple, it was observed in 16J that for Euler's equation, the classical Schwarz algorithm 
with characteristic information exchange at the interfaces is convergent, even without 
overlap. To relate systems of partial differential equations to scalar ones, the algebraic 
tool of the Smith factorization [IT] has proved to be useful, see [IT] . 
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